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Введение. В рамках данного исследования рассмотрены 
прецизионные элементы целевой аппаратуры и чувстви- 
тельные элементы системы стабилизации и ориентации 
перспективных космических аппаратов. Обоснованы и 
разработаны метод и программно-алгоритмическое обес- 
печение моделирования динамических характеристик 
указанных элементов. При этом учтены результаты обра- 
ботки данных по экспериментальным исследованиям ак- 
тивных и пассивных осцилляторов. 

Материалы и методы. Показано, как метод обезвешива- 
ния позволяет воссоздать условия, максимально соответ- 
ствующие реальной эксплуатации перспективных косми- 
ческих аппаратов, прецизионных элементов конструкции, 
целевой аппаратуры и их приводов. Представлены схемы 
соответствующих экспериментальных установок. Исполь- 
зованы методы математического моделирования, методы 
механики и динамики машин. Рассчитаны основные пара- 
метры динамики предлагаемой конструкции — определя- 
ющие с точки зрения реализации целевых функций кос- 
мического аппарата. Сформированы рациональные вари- 
анты компоновки и примерных циклограмм функциони- 
рования перспективных космических аппаратов с целью 
снижения микровозмущений от приводных устройств с 
вращающимися массами. 

Результаты исследования. Обоснован и разработан метод 
моделирования динамических характеристик приводов 
перспективных космических аппаратов с учетом функци- 
онирования штатных осцилляторов. Представлен ком- 
плекс методик решения проблем идентификации динами- 
ческих параметров математической модели перспектив- 
ных космических аппаратов с учетом результатов обра- 
ботки данных, получаемых при экспериментальной отра- 
ботке активных и пассивных осцилляторов. Отмечены два 
вида вибрации от двигателей маховиков. Первый: по ко- 
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мандам системы управления ориентацией и стабилизаци- 
ей. Второй: обусловленный остаточным дисбалансом, от 
измерителя солнечной постоянной. Показано, каким обра- 
зом эти вибрации влияют на динамические характеристи- 
ки посадочных мест гироскопического измерителя векто- 
ра угловой скорости и многозонального сканирующего 
устройства гидрометеорологического обеспечения косми- 
ческого аппарата. Полученные данные предназначены для 
решения задач обеспечения динамической точности пер- 
спективных космических аппаратов. 

Обсуждение и заключения. Предложен метод моделирова- 
ния динамических характеристик перспективных косми- 
ческих аппаратов при функционировании в режиме пре- 
цизионной ориентации. Решение основано на представ- 
ленных В статье результатах теоретико- 
экспериментальных исследований и учитывает работу 
штатных осцилляторов. Реализация указанного метода 
доведена до программно-алгоритмического обеспечения 
оценки динамических характеристик штатных осциллято- 
ров перспективного космического аппарата. Обоснованы 
рекомендации по снижению влияния активных осцилля- 
торов. Выбраны исходные данные для определения дина- 
мики перспективных космических аппаратов с точки зре- 
ния выполнения их целевых функций. Предложены ком- 
поновка и примерные циклограммы функционирования 
перспективного космического аппарата с целью выявле- 
ния приводных устройств с вращающимися массами как 
источников микровозмущений. 


Ключевые слова: амплитуда колебаний, декремент зату- 
хания колебаний, динамическая точность, космический 
аппарат, математическое моделирование, метод, микро- 
возмущение, перемещение, прецизионная стабилизация, 
привод, программно-алгоритмическое обеспечение, ско- 
рость, ускорение, частота колебаний. 
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Введение. На современном этапе развития авиакосмической отрасли сохраняют актуальность практи- 


ческие задачи обеспечения динамической точности угловых и линейных перемещений прецизионных элемен- 
тов конструкции (ПЭК) и целевой аппаратуры (ЦА) перспективных космических аппаратов (КА) при действии 
внутренних источников возмущений. Такими источниками считаются устройства, приборы и приводы косми- 
ческого аппарата и научной аппаратуры. Движущиеся и вращающиеся массы указанного оборудования в ре- 
жимах прецизионного наведения формируют вибрационные возмущения [1,2]. В ходе теоретико- 
экспериментальных исследований решен комплекс перечисленных ниже частных задач исследования [3, 4]. 

1) Анализ требований, предъявляемых к современным космическим аппаратам с прецизионной стаби- 
лизацией, в зависимости от их назначения и устанавливаемой на них научной аппаратуры. Цель — унификация 
этих требований. 

2) Классификация и анализ основных источников внутренних возмущений. Цель — определение 
наиболее виброактивных и выявление возможности ослабления и (или) исключения их воздействия. 

3) Обоснование и разработка экспериментальных методов и средств определения микровозмущений, а 
также методов их математического моделирования. 

4) Проведение экспериментов по определению диссипативных свойств конструкций космических ап- 
паратов при малых уровнях перемещений (порядка 0,5 мкм). 

5) Обоснование и разработка достаточно точных математических моделей космических аппаратов и 
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научной аппаратуры для анализа динамической точности как в низкочастотной, так и в среднечастотной обла- 
стях возмущений. 

6) Обоснование и разработка методов математического моделирования для исследования динамиче- 
ской точности с учетом экспериментально определяемых параметров возмущений. 

7) Анализ результатов исследования динамической точности и общих требований к ее параметрам для 
перспективного космического аппарата, включая ЦА. Обоснование и разработка на основе этого анализа част- 
ных требований по виброактивности основных источников внутренних возмущений. 

Материалы и методы 

Постановка задачи исследования. Предложены схемы экспериментальных установок, реализация ко- 
торых с применением метода обезвешивания позволяет воссоздать условия, максимально соответствующие 
реальной эксплуатации перспективных КА, ПЭК, ЦА и их приводов [3—5]. Математическое моделирование и 
эксперименты позволили обосновать и разработать предложения по коррекции динамических характеристик 
ПЭК и ЦА перспективных КА при действии внутренних источников возмущений. 

Методы решения задачи исследования. Для решения задач исследования использованы методы ма- 
тематического моделирования, методы механики и динамики машин. Предложены: метод моделирования ди- 
намических характеристик перспективных КА при функционировании в режиме прецизионной ориентации с 
учетом работы штатных осцилляторов и метод обработки результатов экспериментальных исследований мик- 
ровозмущений. 

Описание схем экспериментальных установок. Для решения указанных задач реализованы схемы 
проведения экспериментальных исследований (рис. 1). 
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Рис. 1. Схемы проведения эксперимента на силовом испытательном стенде с нагрузкой на выходном валу модуля привода: 


горизонтальная ось вращения (4); вертикальная ось вращения (6): нагрузка на выходном валу ./= 4,75 кгм* (в) 


Экспериментальные исследования по выбранным схемам позволили оценить влияние гравитации и си- 
стемы обезвешивания на полезный сигнал [5-8]. 

Результаты исследования 

Метод моделирования динамических характеристик приводов перспективных космических ап- 
паратов с учегом функционирования штатных осцилляторов. Для обеспечения динамической точности 
принят представленный на рис. 2 порядок проведения исследований [8—11]. 


Испытания источников возмущений 
с целью определения их силового 


Создание мат. модели прохождения 
1 микроколебаний ОТ ИСТОЧНИКОВ 


х х и моментного воздействия 
возмущении к местам прецизионнои ЦА 


на конструкцию КА 


Расчет динамической реакции ЦА 
2 на воздействие микровозмущений 





Оценка качества работы ЦА Выработка требований по допустимым 
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З воздействием источников возмущений крепления приборов 
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5 до достижения необходимых 


характеристик 





Рис. 2. Порядок проведения исследований, принятый при отработке динамической точности 


Кроме того, экспериментально исследованы микровозмущения физических натурных моделей перспек- 
тивных КА. При анализе информации, получаемой при экспериментальной отработке активных и пассивных 
осцилляторов, возникают проблемы идентификации динамических параметров математической модели изде- 
лия. В данной статье представлены пути и методы решения этих проблем. 

Метод обработки результатов экспериментальных исследований микровозмущений перспектив- 
ных КА. Метод определения основных параметров динамической схемы основывается на предположении, что 
полученная выборка свободного движения представляет собой общее решение системы линейных дифференци- 
альных уравнений, то есть имеет вид [3, 4, 7, 8]: 
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И 
—01. 
у =» де" зт(2/и+Ф,), (1) 
= 
где и — общее число искомых составляющих в выборке; А; — амплитуда вклада /-го тона; О ; —_ декремент 


затухания /-го тона; // — частота /-го тона колебаний; ф, — фаза /-го тона колебаний; { — время. 


Необходимо максимально приблизить эту зависимость к колебательной части сигнала с прибора, полу- 
ченного по каналам телеметрии [7, 8]. Для фиксации идентифицируемого сигнала нужно провести следующие 
предварительные операции [3, 4, 7, 8, 10, 11]: 

— выделить участок телеметрической информации по угловым скоростям или углам, где наблюдаются 
колебания, определяемые влиянием жидкого наполнения; 

— удалить из полученного сигнала низкочастотные и постоянные составляющие, обусловленные угловым 
движением космического аппарата относительно его центра масс (производится путем фильтрации 
низкочастотной составляющей и (или) удалением полиномиального тренда, соответствующего твердотельным 
формам движения изделия); 

— привести сигнал к постоянному интервалу регистрации; 

— определить частотные диапазоны основных гармонических составляющих сигнала на основе использования 
алгоритма быстрого преобразования Фурье; 

— сместить сигнал, чтобы время первого замера в исследуемой части сигнала соответствовало началу 
временной оси (это необходимо для правильной оценки диссипативных и фазовых характеристик). 

Преобразованный таким способом сигнал используется в качестве таблично заданной функции (х;) при 
его аппроксимации зависимостью для определения параметров процесса. При подборе параметров 


используется метод наименьших квадратов. Параметры А, 6 РЕ Ди ф ; определяются путем минимизации [3, 4, 


7, 8]: 


т И 
—067хЕ <: 2 : 
У -> де зш (2/7; +ф,;))” > ша, (2) 
1-1 75 
где х; — значения сигнала, полученного при обработке результатов испытаний; # — моменты времени, 
соответствующие замерам х; т — количество замеров; 4; — амплитуда вклада /-го тона; 0; — декремент 


затухания /-го тона; /— частота /-го тона колебаний; ф, — фаза /-го тона колебаний. 


Результаты моделирования и экспериментальные исследования учитывают следующие основные 

факторы: 

— количество гармонических составляющих, проявляющихся в исследуемой выборке; 

— близость расположения частот этих составляющих; 

— длина выборки (определяет и позволяет достоверно оценить необходимое число периодов движения по 
самой низкой частоте); 

— шаг записи сигнала (определяет необходимое число точек на периоде движения самой высокой частоты); 

— разброс амплитуд отдельных движений составляющих выборки; 

— наличие в выборке ненаблюдаемых входных воздействий. 

Разумеется, все перечисленные факторы взаимозависимы, а их влияние на точность и достоверность 
результата определяется их сочетанием. Поэтому исследование проводится для каждого фактора при вариации 
остальных. 

Созданный метод апробирован при работе осцилляторов, в результате которой силовые факторы воз- 
действуют на места крепления. При определении амплитудно-частотных характеристик (АЧХ) сил и моментов 
(работа модуля на штатной частоте импульсов шагового двигателя (ШД) 130 Гц) с учетом кинематической цепи 
конструкции применялась следующая последовательность разработки модуля: 

— выбор из массива данных наиболее характерных результатов измерений (критерии выбора — максимальные 
амплитуды измеренных силовых параметров и максимальный набор пиков в АЧХ силовых параметров); 

— анализ результатов измерений внешнего фона и колебаний систем крепления и обезвешивания в частотной 
области с целью выявления частот шумовых сигналов; 

— составление списка диапазонов шумовых частот для каждого силового параметра; 

—щ фильтрация из исходных процессов изменения силовых параметров в выбранных результатах измерения 
шумовых диапазонов частот (используется специальный пакет прикладных программ); 

—щ сохранение результатов измерения отфильтрованных процессов и анализ их спектров с целью выявления 
основных осцилляторов кинематической цепи модуля. 
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Отдельные этапы разработки модуля описаны в [3, 4, 7, 8, 10, 11]. 

По представленным на рис. | схемам были проведены измерения сил и моментов, приходящих на места 
крепления модуля, для трех вариантов нагрузки на выходной вал со значениями момента инерции: 
1] = 1,53 кгм"; 2) Л, = 6,1 кгм”; 3) Л = 7,6 кгм°. 

Цель экспериментов — определить, как максимальные значения крутящего момента относительно оси 
вращения выходного вала зависят от величины момента инерции нагрузки на валу, а затем составить прогноз 
по величине крутящего момента для штатной нагрузки. При определении этой зависимости и крутильной жест- 
кости на выходном валу модуля к анализу принимались измерения момента М. относительно оси вращения вы- 
ходного вала. 

Крутильная жесткость определялась по формуле: 

с= 7.4? Г Л (3) 
где Л. — момент инерции нагрузки на выходном валу, кгм”; }; — низшая частота крутильных колебаний относи- 
тельно оси вращения выходного вала, Гц; { — номер варианта нагрузки на выходной вал модуля. 

Из полезного сигнала момента М. (фильтрации) выделили процесс, в котором изменение данного мо- 
мента происходит на низшей собственной частоте крутильных колебаний. Это позволило определить декремен- 
ты колебаний. 

Декременты крутильных колебаний выходного вала на низшей собственной частоте с различными ва- 
риантами нагрузки определялись по формуле: 


6 = С, я р (4) 


где 4„ — амплитуда ускорения на т-м цикле колебаний; А» + „ — амплитуда ускорения на (т - п)-м цикле коле- 
баний; ип — число принятых к анализу циклов колебаний. 

При определении жесткости выходного вала и декрементов крутильных колебаний на низшей соб- 
ственной частоте применялась описанная ниже последовательность действий. 

1) Выбор трех наиболее характерных процессов изменения момента М. для трех вариантов нагрузки на 
выходном валу модуля. 

2) Определение низшей собственной частоты крутильных колебаний. В этом случае по максимальной 
амплитуде в АЧХ момента М. определяется частотный диапазон, в котором предположительно находится зна- 
чение низшей собственной частоты крутильных колебаний. С помощью узкополосной фильтрации анализиру- 
ются узкие диапазоны частот в выбранной области (критерием низшей собственной частоты крутильных коле- 
баний среди анализируемых отфильтрованных процессов является монотонное убывание амплитуды во време- 
ни и отсутствие биений). 

3) По значениям низшей собственной частоты крутильных колебаний и момента инерции нагрузки с 
использованием формулы (3) определяется жесткость выходного вала. 

4) Из полученного временного процесса изменения амплитуд крутильных колебаний на низшей соб- 
ственной частоте определяются значения Аи, Аш + „, п, а по формуле (2) — декремент этих колебаний. 

Особенности обработки результатов исследований воздействия на места крепления со стороны 
силовых факторов, возникающих при работе осцилляторов. Автоматически выполняется преобразование 
формата первичных данных измерений в формат, удобный для графического представления и анализа инфор- 
мации. Применение специальной программы позволяет значительно сократить время первичной обработки ре- 
зультатов измерений. Для исключения из исходного процесса шумовых сигналов также используется пакет 
специальных прикладных программ. 

Для каждого из измеренных силовых параметров (К, Р’, Р., М., М,, М.) характерны не только общие 
для всех них шумовые сигналы (например, сетевая электрическая помеха на частоте —-50 Гц или стабильное 
внешнее возмущение грунта в диапазоне частот от 16 до 17,5 Гц), но и собственные, обусловленные влиянием 
систем обезвешивания и крепления привода. Для каждого из силовых параметров требуется исключить в сред- 
нем двенадцать частотных диапазонов. Соответственно, в каждой выбранной для анализа совокупности резуль- 
татов измерений требуется исключение разных наборов шумовых сигналов для каждого из силовых парамет- 
ров, если решается задача определения динамических характеристик модуля с учетом кинематической цепи 
(при этом необходим анализ измеренных силовых параметров в широком диапазоне частот от 0 до 200 Гц). 
Критериями выбора результатов измерений силовых факторов являются: максимальные амплитуды силовых 
параметров; максимальный набор пиков в амплитудно-частотных характеристиках силовых параметров. 

При обработке результатов измерений в задаче определения жесткости модуля на выходном валу и де- 
кремента колебаний на низшей собственной частоте необходимо из исходного сигнала (фильтрации) выделить 
процесс, в котором изменение силового параметра происходит на низшей собственной частоте крутильных ко- 
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лебаний. Эта процедура успешно реализуется в специальном пакете прикладных программ [3-4]. 

Методика выбора характерных режимов испытаний и соответствующих им результатов измере- 
ний. Основные схемы испытаний и режимы работы (частота импульсов шагового двигателя (ШД) 130 Гц) мо- 
дуля были выбраны с учетом максимального приближения к условиям функционирования КА на орбите. Кри- 
териями выбора наиболее характерных режимов для решения поставленной задачи являются: 

— отсутствие нагрузки на выходном валу; 
— вертикальная ось вращения выходного вала; 
— частота опроса данных измерений не менее 500 Гц (исследуемый диапазон частот — от 0 до 200 Гц). 

Критериями выбора результатов измерений для определения динамических характеристик модуля с 
учетом кинематической цепи являются: 

— максимальные амплитуды измеренных силовых параметров; 
— максимальный набор пиков в АЧХ силовых параметров. 

Методика выбора характерных режимов и результатов измерений для определения жесткости и 
декремента крутильных колебаний выходного вала модуля. Критериями выбора характерных режимов для 
решения поставленной задачи являются: 

— вертикальная ось вращения выходного вала (см. рис. 1,6); 
— три варианта нагрузки на выходном валу (величина нагрузки не менее 0,5 кгм”); 
— частота опроса данных измерений не менее 50 Гц (исследуемый диапазон частот — от 0 до 10 Гц). 

Критерием выбора файлов являются максимальные амплитуды измеренного момента М... 

Методика определения динамических характеристик модуля с учетом кинематической цепи. По- 
сле обработки исходного измеренного сигнала формируется полезный сигнал, отражающий динамику звеньев 
кинематической цепи модуля при работе ШД на частоте управляющих импульсов 130 Гц. Анализ АЧХ каждого 
из силовых параметров (К,, К, ЁР., М,, М, М.) позволил выявить основные диапазоны частот, на которых имели 
место повышенные амплитуды этих параметров. В результате проведенных исследований определены времен- 
ные процессы и спектры для сил Р,, РЁ’, РЁ, и моментов М, М,, М... На рис. 3 представлены примеры временных 
процессов и спектров для силы РЁ’, и момента М... 
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Рис. 3. Временной процесс (а) и спектр (6) силы ЁЕ,; временной процесс (в) и спектр (г) момента М, на частоте ШД 130 Гц 


Результаты определения динамических характеристик модуля на низшей собственной частоте 


крутильных колебаний выходного вала. Результаты расчетов приведены в табл. 1. 
Таблица 1 


Частотные характеристики по параметрам сил Р,, К’, Р, и моментов М,, М, М, при работе модуля на частоте 
импульсов ШД 130 Гц 


№ п/п | Параметр нагрузки | Частота /, Гц (в скобках указано осредненное значение амплитуды нагрузки) 


0,5-0,9 (0,014); 5=6,1 (0,08); 8,2 (0,04); 12,3 (0,03); 26/2 (0,04); 45 (0,022); 
91(0,02); 105 (0,014); 130 (0,39); 170 (0,01); 205 (0,04) 
56,1 (0,014): 8,2 (0,02):12,3 (0,02); 26,2 (0,044); 45 (0,014); 91 (0,017): 130 
(0,47); 170 (0,02); 205(0,04) 
0,5-0,9 (0,07); 26,2 (0,04): 45 (0,015); 91 (0,04); 130 (0,1): 170 (0,02) 


0,5-0,9 (0,01); 8,2 (0,0025); 12,3 (0,03): 26,2 (0,025): 91 (0,0026): 105 (0,003); 130 





: (0,3); 170(0,001); 205 (0,0045) 

Е 56,1 (0,016); 8,1 (0,0062); 26,2 (0,016); 45 (0,004); 91 (0,008); 105 (0,005); 

Е 130(0,075); 170 (0,01); 205 (0,0055) 

ги 0,5-0,9 (0,0025); 5=6,1 (0,005); 18,5 (0,004); 26,2 (0,011); 45 (0,006); 105 

- (0,0018); 130 (0,3); 170 (0,015) 

_ 

. После обработки исходных измеренных сигналов по предложенной методике получены данные для 

Е определения жесткости выходного вала по формуле и декрементов крутильных колебаний. Анализ крутильных 
колебаний выходного вала модуля проводился для моментов инерции: 1,53 кг-м”; 6,1 кг.м”; 7,6 кг-м°. При ис- 
следованиях крутильных колебаний выходного вала модуля сравнивались две ситуации: 

324 — ШД включен и кинематическая цепь в работе; 


— ШД выключен. 
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В первом случае низшая собственная частота крутильных колебаний ниже, чем во втором. Результаты 
анализа приведены в табл. 2. 
Таблица 2 
Результаты определения жесткости выходного вала модуля с нагрузкой на низшей собственной частоте 


Момент Низшая собственная частота системы, Гц Жесткость выходного вала, Н`м/град 


инерции 
нагрузки, кг.м 


[6 


Я О Е О ПЕ 
ПР ОИ ПОТ ОИ ПОТ О ОО ОО ПО ООИ 


Полученные осредненные величины амплитуд моментов №2 и декрементов колебаний для низших соб- 





ственных частот системы приведены в табл. 3. 
Таблица 3 
Декременты крутильных колебаний и осредненные амплитуды момента относительно оси вращения выходного вала 


Примеры временного процесса момента М2 и его спектра для случая нагрузки на выходном валу моду- 

ля .Л› = 7,6 кгм” представлены на рис. 4. 
М2, Нм 
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Рис. 4. Графики временного процесса (а) и спектра (6) момента М2 для случая нагрузки на выходном 
валу модуля .Л› = 7,6 кгм” 
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Обсуждение и заключения. Исследования позволили определить динамику приводов перспективного 
КА и его составных частей как источников внутренних возмущений. Анализ результатов исследований показал, 
что для снижения влияния активных осцилляторов на динамическую точность перспективного КА в режиме 
прецизионной ориентации необходимо: 
— провести экспериментальное подтверждение исходных данных по бортовым источникам возмущений; 
—щ провести экспериментальное исследование диссипативных характеристик и характеристик жесткости кон- 
струкции перспективного КА; 
— разработать критерии оценки влияния воздействий от бортовых источников возмущений на целевое функ- 
ционирование перспективного КА в диапазоне частот до 100 Гц; 
— уточнить среднечастотную динамическую модель перспективного КА по результатам экспериментальных 
исследований; 
—Й провести анализ чувствительности бортовых приборов перспективного КА к вибрациям и разработать тре- 
бования к виброактивности источников микровозмущений; 
—щ провести определение виброактивности источников вибровозмущений: двигателя-маховика (ДМ), привода 
зеркала МСУ-ГС, измерителя солнечной постоянной ИСП-2М. 

Результаты проведенных теоретико-экспериментальных исследований позволили решить перечислен- 
ные ниже задачи. 

1) Анализ источников микровозмущений: 
— анализ состава, режимов работы и характеристик бортовых источников микровозмущений; 
— оценка вибровозмущений конструкции КА, генерируемых электромеханическими исполнительными орга- 
нами (ЭМИО); 
— оценка моментов, создаваемых ДМ по командам системы управления ориентацией и стабилизацией (СУОС); 
— оценка вибровозмущений конструкции КА, генерируемых многозональным сканирующим устройством гид- 
рометеорологического обеспечения (МСУ-ГС); 
— оценка вибровозмущений конструкции КА, генерируемых приводом солнечной батареи (СБ); 
— оценка вибровозмущений конструкции КА, генерируемых приводом остронаправленной антенны (ОНА); 
—щ оценка вибровозмущений конструкции КА, генерируемых измерителем солнечной постоянной величи- 
ны (ИСП-2М). 

2) Оценка динамических характеристик посадочных мест гироскопического измерителя вектора угло- 
вой скорости (ГИВУС) перспективного КА при вибрационных воздействиях от ИСП-2М, а также от: 
— МСУ-ГС, 
— ДМ (вибрации обусловлены командами СУОС и остаточным дисбалансом). 

3) Оценка динамических характеристик посадочного места МСУ-ГС перспективного КА при возмуще- 
ниях от ДМ по командам СУОС, ДМ, обусловленных остаточным дисбалансом, ИСП-2М. 
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